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1  Optimization  of  Cloaks  and  Anti-Reflection  Coatings 


We  have  simulated  a  2D  multilayer  cloak  wrapped  around  a  perfectly  magnetic  conducting  cylinder 
with  an  incident  TM  plane  wave  as  shown  in  Fig.  1.  Each  cylindrical  layer  of  the  cloak  in  Fig.  1(a) 
has  piecewise  constant  material  parameters  er,  e §  and  jiz.  The  material  values  were  obtained  using 
our  optimization  procedure  and  by  sampling  the  continuous  material  parameters  of  Pendry’s  original 
2D  cylindrical  cloak  at  the  middle  of  each  layer.  The  optimization  was  achieved  using  a  nonlinear 
constrained  optimization  procedure  that  was  coupled  with  a  forward  finite  element  solver,  with  a 
cost  function  c(er,e0,^_)  defined  as  an  azimuthal  integration  over  the  radar  cross  section  (RCS)  of 
the  cloaked  cylinder, 


)  =  lim 
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2tt  2 TTy/r  Es(r,0;er,e^,^) 
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where  er,  ,  and  /jy  are  vectors  containing  the  respective  constitutive  parameters  of  each  of  the 
cloak  layers,  r  and  (j)  denote  the  radial  and  azimuthal  coordinates  respectively,  Es(r,  0;  er,  e^,  fi  )  is  the 
scattered  field  of  a  cloak  defined  by  (er,f0,  [ij  at  the  coordinate  (r,  </>),  and  E*  is  the  incident  field. 
Fig.  1(b)  shows  the  RCS  of  the  cloaked  cylinder  both  before  and  after  optimization  for  different 
number  of  layers  in  the  cloak  discretization.  It  is  evident  from  this  figure  that  the  optimization 
procedure  improves  the  performance  of  the  cloak  even  when  it  is  discretized  into  a  small  number 
of  layers,  but  as  the  number  of  layers  is  increased,  the  correspondingly  larger  number  of  degrees 
of  freedom  in  the  optimization  helps  in  achieving  substantial  improvement.  Fig.  1(c)  gives  the 
simulated  magnetic  scattered  field  from  a  cloaked  cylinder  comprising  8  layers  both  before  and  after 
optimization.  This  figure  clearly  shows  that  the  optimized  cloak  induces  reduced  scattering  even  in 
the  near  field. 

Using  the  complex  constraint  for  the  minimization,  another  cloak  optimization  has  been  studied. 
These  constraints  are  related  with  physical  fabrication  limits  and  lead  to  cloak  designs  that  can 
be  realized  using  current  fabrication  tools  and  materials.  An  example  for  such  a  constraint  is  the 
dynamic  range  of  the  anisotropic  material  properties  that  could  be  obtained  using  layered  metal- 
insulator  stacks,  a  structure  we  have  been  studying  recently  [1],  An  initial  estimate  for  this  dynamic 
range  has  been  obtained,  assuming  effective  medium  conditions  and  a  large  cloak  radius  of  curvature, 
from  the  dielectric  constants  of  a  1-D  metal-insulator  stack.  The  frequency-dependent  anisotropic 
dielectric  constants  for  an  alternating  stack  of  silver  and  magnesium  oxide  (MgO)  for  different  metal 
duty  cycles  D  are  shown  in  Fig.  2(a).  From  Fig.  2(a),  the  dispersion  will  be  smallest  (and  thus, 
the  cloaking  performance  will  be  relatively  wideband)  if  operation  is  beyond  resonance,  i.e.,  in  the 
longer  wavelength  range.  This  regime  gives  a  dielectric  constant  parallel  to  the  stack  interface  (eii) 
which  is  smaller  than  that  perpendicular  to  the  interface  (ej_).  Keeping  the  cloak  non- magnetic,  good 
performance  can  be  achieved  with  Ag  and  MgO  layers  stacked  azimuthally  (see  Fig.  2(b)),  thereby 
achieving  a  radial  dielectric  constant  that  is  smaller  than  the  azimuthal  dielectric  constant.  The  next 
step  toward  realizing  an  optical  cloaking  device  will  be  to  use  the  bounds  on  the  dynamic  range  of  the 
anisotropic  material  properties  of  the  cloak  envisioned  in  Fig.  2(b)  as  constraints  for  the  optimization 
procedure  aimed  at  the  design  of  an  efficient  cloak.  To  exemplify  the  cloak  performance  using  the 
suggested  optimization  procedure,  even  when  the  initial  material  parameters  stem  from  the  reduced 
nonmagnetic  cloak  parameter  set  and  the  optimization  is  constrained  to  keep  the  cloak  nonmagnetic, 
Fig.  2(c)  shows  the  RCS  of  the  nonmagnetic  cloaked  cylinder  both  before  and  after  optimization  for 
8  layers  in  the  cloak  discretization.  The  dramatic  improvement  in  the  cloaking  performance  when 
optimizing  a  cloak  based  on  the  nonmagnetic  reduced  material  properties  [2]  encourages  us  to  believe 
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that  the  implementation  of  the  bounds  of  the  dynamic  range  of  metal-insulator  stacks  as  constraints 
for  the  optimization  procedure  will  yield  effective  and  realistic  cloaking. 


(a)  (b)  (c) 

Figure  1:  (a)  A  two-dimensional  N  layer  cloak  wrapped  around  a  perfectly  magnetic  conducting 
cylinder  with  an  incident  TM  plane  wave  impinging  from  the  left,  (b)  Simulated  RCS  of  a  cloaked 
cylinder  with  inner  radius  of  R\  =  700  nm  and  outer  radius  of  R2  =  1400  nrn,  illuminated  at 
A  =  400  nm.  The  cloak  is  discretized  into  4,  6  and  8  layers,  each  shown  in  a  different  panel  in  which 
the  continuous  blue  line  corresponds  to  the  result  by  discretizing  Pendry’s  original  2D  cloak  [3] 
and  the  doted  red  line  corresponds  to  cloak  parameters  obtained  using  optimization,  (c)  Simulated 
scattered  magnetic  field  for  the  8  layer  cloak  with  the  same  geometry  and  incidence  as  in  (b),  before 
(top)  and  after  (bottom)  the  optimization  process.  The  total  field  inside  the  cloaked  region  is  zero 
because  the  scattered  field  (shown)  cancels  the  incident  held. 


(a)  (b)  (c) 

Figure  2:  (a)  Calculated  real  dielectric  constants  en  and  ej_,  parallel  and  perpendicular  to  the  surface 
of  the  effective  slab,  for  a  Ag/MgO  stack  with  Ag  volume  fraction  D,  using  published  data,  (b) 
Conceptual  illustration  of  a  cloak  built  out  of  silver  and  magnesium  oxide  stacks,  (c)  Simulated  RCS 
(dB)  of  a  700  nm  thick  nonmagnetic  cloak  wrapped  around  a  PMC  cylinder  of  1400  nm  in  diameter 
that  is  illuminated  at  A  =  400  nm.  The  cloak  is  discretized  into  8  layers.  The  continuous  blue  curve 
corresponds  to  the  RCS  of  only  the  PMC  cylinder  without  a  cloak,  the  dashed  red  curve  corresponds 
to  the  RCS  of  a  PMC  cylinder  covered  by  a  cloak  with  parameters  obtained  by  discretizing  the 
reduced  nonmagnetic  2D  cloak  suggested  in  [2]  and  the  doted  green  curve  corresponds  to  a  PMC 
cylinder  covered  by  a  nonmagnetic  cloak  with  optimized  parameters. 
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(a) 


(b) 


(c)  (d) 

Figure  3:  Power  transmission  (T),  reflectance  (R),  and  absorption  (A)  for  a  2  /i.m  thick  planar  slab 
composed  of  QD/Ta  mixture  as  a  function  of  the  incident  wavelength  A  for  varying  excited  QD 
fraction.  The  percentage  of  excited  QDs  were  chosen  to  be  (a)  0%  (b)  25%  (c)  37.5%  and  (d)  55%. 
The  fill  fraction  of  tantalum  (Ta)  and  CdSe  QDs  in  an  SiC>2  background  is  0.16  and  0.04,  respectively. 
Note  the  opening  of  the  transmission  window  in  the  slab  with  the  increase  in  the  percentage  of  the 
excited  QDs  in  the  mixture. 


2  Active  Quantum  Dot  Optical  Metamaterials 

Realizing  an  active  nretamaterial  composed  of  a  semiconductor  quantum  dot  (QD)  and  metal  nanopar¬ 
ticle  mixture  provides  important  frequency-selective  gain  and  loss  properties,  depending  on  size  and 
composition,  and  with  optical  (or  ultimately  electrical)  pumping,  the  metamaterial  response  can  be 
controlled.  For  example,  it  is  in  principle  possible  to  achieve  a  lossless  metal,  and  this  would  revolu¬ 
tionize  plasmonics  and  open  new  domains  for  ultrafast  signal  processing  and  sensing.  The  material 
is  also  key  to  nanophotonic  sources  and  detectors,  because  of  compatibility  of  fabrication  with  plas- 
monic  materials.  Particularly  useful  is  the  control  over  the  loss  and  gain  with  the  narrow  exciton 
level  in  the  QD,  making  possible  dispersion  engineering  and  a  lossless  local  metallic  response.  We  are 
pursuing  an  experimental  effort  to  demonstrate  QD  metamaterials  with  important  characteristics  in 
applications. 

Our  preliminary  studies  show  that  using  the  properties  of  an  active  QD  mixture,  it  is  possible 
to  achieve  a  material  response  with  a  narrowband  transmission  window  while  maintaining  low  re¬ 
flectance.  Figure  3  shows  the  calculated  power  transmission  (T),  reflectance  (R),  and  absorption  (A) 
of  an  example  QD/metal  mixture  slab,  placed  in  vacuum,  for  different  fractions  of  excited  QDs  in 
the  mixture.  This  result  uses  tantalum  (Ta)  as  the  metal,  offering  high  broadband  absorption,  and 
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Figure  4:  Shift  in  the  central  wavelength  of  the  power  transmission  window  with  the  change  in  the 
diameter,  d,  of  the  QDs  used  in  the  gain  mixture  in  the  2  /xm  thick  planar  slab.  The  fill  fraction  of 
tantalum  (Ta)  and  CdSe  QDs  in  SiC>2  background  was  0.16  and  0.04,  respectively.  The  percentage 
of  excited  QDs  was  chosen  to  be  55%  for  each  case. 

CdSe  QDs  as  the  gain  material  to  open  up  the  transmission  window.  The  fill  fraction  of  Ta  and  CdSe 
QDs  in  an  SiC>2  background  for  all  the  cases  presented  is  0.16  and  0.04,  respectively.  The  QDs  have 
a  diameter  of  1.9  nrn  and  the  dielectric  constant  was  obtained  using  our  simplified  QD  permittivity 
model  [4],  The  effective  permittivity  for  the  mixture  was  then  obtained  using  the  Maxwell  Garnett 
mixing  formula  with  three  species  (tantalum,  gain  QDs  and  loss  QDs)  in  the  Si02  background  [4]. 
The  power  transmission  and  reflection  for  the  planar  isotropic  slab  is  calculated  using  the  standard 
transmission  and  reflection  formulas  with  a  normally  incident  plane  wave  [5].  From  Fig.  3,  a  trans¬ 
mission  window  opens  as  the  pump  energy  is  increased  or  with  an  increase  in  the  fraction  of  the 
excited  QDs  in  the  mixture.  In  this  example  of  a  2  //m  thick  slab,  the  QDs  provide  the  gain  window 
and  the  Ta  nanoparticles,  which  constitute  16.0  %  of  the  mixture,  absorb  most  of  the  energy  at  other 
frequencies.  The  overall  effect  of  the  mixture  allows  us  to  open  a  narrow  transmission  window  by 
tuning  the  excited  QD  concentration  using  different  pump  energies.  Figure  4  shows  the  shift  in  the 
central  wavelength  of  the  transmission  window  with  the  change  in  the  diameter,  d,  of  the  QDs  used 
in  the  gain  mixture  of  the  2  /xm  thick  planar  slab.  The  fill  fraction  of  tantalum  and  CdSe  QDs  in  the 
SiC>2  background  was  again  0.16  and  0.04,  respectively,  and  55%  of  QDs  were  assumed  to  be  excited. 
The  center  of  the  transmission  window  shifts  to  longer  wavelength  with  an  increase  in  the  diameter 
of  the  QDs.  A  proper  choice  of  the  QD  diameter  and  the  metal  nanoparticle  can  theoretically  allow 
us  to  open  up  a  transmission  window  at  any  specific  wavelength. 

Another  application  for  active  QD  metamaterials  is  as  the  gain  medium  in  a  plasmonic  amplifier. 
Surface  plasmon  waves  at  metal  surfaces  allow  light  to  be  tightly  confined,  and  various  waveguide 
configurations  provide  great  potential  for  future  optical  devices.  With  active  optical  metamaterials 
having  low  loss  or  even  overall  gain,  the  propagation  length  of  the  surface  plasmon  waves  can  be 
significantly  enhanced.  Also,  with  an  active  metamaterial  having  overall  gain  inside  a  plasmonic 
cavity,  a  nano-laser  can  be  achieved.  We  have  studied  the  use  of  a  quantum  dot  (QD)  and  dielectric 
mixtures  in  a  metal  slot  cavity  to  create  a  plasmonic  amplifier.  We  have  previously  reported  that  a 
small  resonant  cavity,  with  the  gap  size  much  smaller  than  wavelength,  is  capable  of  confining  the  light 
and  providing  highly  localized  resonant  field  enhancement  for  specific  nanoparticle  structures  [6, 7] . 
The  cavity  resonant  length  can  be  predicted  using  the  wavelength  of  the  dominant  propagating 
waveguide  mode.  The  cavity  quality  factor  ( Q )  is  dictated  by  the  geometry  and  the  complex  dielectric 
constant  of  the  metal.  In  designing  and  simulating  active  plasmonic  devices,  the  material  properties, 
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Figure  5:  (a)  Schematic  view  of  a  MIM  waveguide  cavity  amplifier  with  a  CdSe/Si02  QD  mixture, 
(b)  FEM  simulation  of  the  total  power  flow  magnitude  [W/m2]  for  red  light  (A=  633  nm).  The 
absorption  coefficient  for  these  QDs  in  a  Si02  background  was  found  to  be  1.1  x  107  m_1  at  the 
pump  wavelength  of  532  nm  [8]  with  a  QD  fill  fraction  of  0.3.  Using  this  and  a  pump  volume  of 
1.57  x  10-1'  m3  requires  0.17  W  to  pump  60%  of  the  QDs. 


especially  of  the  gain  medium,  are  very  important.  Our  previous  work  on  modeling  the  dielectric 
constant  of  QD  mixtures  [4]  played  a  key  role  in  this  work.  We  used  Maxwell  Garnett  mixing  with 
a  model  for  the  QD  dielectric  constant,  along  with  experimental  QD  data  for  commercially  available 
QDs.  The  QD  amplifier  we  consider  uses  a  (metal-insulator-metal)  MIM  waveguide,  and  a  schematic 
of  the  structure  is  shown  in  Fig.  5(a).  The  cavity  Q  has  been  enhanced  by  placing  two  partially 
reflecting  metal  mirrors  at  each  end  of  the  cavity,  as  shown  in  Fig.  5(a).  An  finite  element  method 
(FEM)  simulation  was  carried  out  to  study  the  amplification  from  such  an  active  cavity,  and  the 
magnitude  of  the  total  power  flow  for  the  amplifier  is  shown  in  Fig.  5(b).  The  fill  fraction  of  QDs 
(having  diameter  2.1  nm)  in  the  Si02  background  is  0.3,  and  the  homogenized  effective  permittivity 
of  the  mixture  is  from  Maxwell  Garnett  mixing  with  two  species  (gain  and  loss  QDs).  A  series  of 
such  FEM  simulations  with  varying  wavelength  A  illustrates  the  amplification  inside  the  cavity  as 
a  function  of  the  incident  signal  wavelength,  and  the  results  are  shown  in  Fig.  6.  The  fraction  of 
excited  QDs  in  the  mixture  is  dependent  on  the  pump  energy,  which  can  be  tuned  to  obtain  different 
homogenized  dielectric  constants  for  the  QD  mixture.  From  Fig.  6,  it  is  clear  that  as  we  increase  the 
excited  QD  concentration,  the  amplification  of  the  signal  in  the  cavity  increases. 

These  examples  illustrate  the  potential  of  QD  metamaterials.  We  have  proposed  an  idealized 
model  for  the  dielectric  constants  of  active  QDs  (see  [4]),  and  now  we  are  pursuing  experiments  to 
better  characterize  the  QDs  by  matching  our  model  parameters  to  the  experimental  data.  This  will 
place  us  in  a  position  to  be  able  to  accurately  design  material  responses  for  several  key  demonstrations 
that  will  open  new  application  domains. 
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Figure  6:  Total  magnitude  of  the  emitted  power  from  an  amplifier  cavity  (shown  in  Fig.  5(a))  as 
a  function  of  the  incident  wavelength  A  for  different  fractions  of  excited  QDs  in  the  mixture.  The 
simulation  parameters  remained  the  same  as  for  Fig.  5(a)),  except  that  the  percentage  of  excited 
QDs  in  the  mixture  is  varied.  The  fraction  of  the  excited  QDs  can  be  tuned  by  tuning  the  pumping 
energy. 


3  Black  Materials  Based  on  Graphene  Stacks 

Relying  on  recent  developments  in  theory  and  measurement  of  the  properties  of  graphene  sheets 
[9-11],  we  have  been  able  to  calculate  the  effective  medium  parameters  for  a  one  dimensional  stack 
of  graphene  sheets  separated  by  an  arbitrary  substance  (see  Fig  7(a))  [12].  As  the  surface  conductivity 
of  a  single  graphene  sheet  is  very  close  to  the  universal  value  of  <7o  =  7re2/2 h  throughout  a  substantial 
portion  of  the  optical  regime  (250  nm-6000  nm),  the  effective  parameters  of  the  stack  vary  linearly 
and  rather  slowly  with  frequency  for  visible  light.  A  few  conditions  must  hold  in  order  for  the 
electromagnetic  homogenization  to  be  valid  [13].  First,  the  distance  between  the  graphene  sheets 
must  be  large  enough  compared  to  the  interatomic  distance  (~  0.3  nm)  to  avoid  mutual  effects 
causing  a  change  in  the  electronic  structure  of  the  sheets.  Two  other  conditions  stem  from  the 
need  for  a  meaningful  homogenization  of  the  stack  bounding  both  the  distance  between  the  sheets 
and  the  conductance  of  each  sheet  to  be  small.  We  show  that  for  visible  and  near-infrared  light 
incident  on  a  stack  immersed  in  a  low  refractive  index  medium,  all  those  conditions  hold  when  the 
distance  between  the  sheets  is  in  the  range  1  nm  <  d  <  Ao/10  nm.  A  low  refractive  index  medium 
is  crucial  for  the  performance  of  such  a  graphene  stack  as  a  dark  material,  and  we  assume  that 
due  to  the  remarkable  mechanical  strength  of  the  sheets,  it  is  possible  to  use  a  sparsely  populated 
matrix  of  small  and  low  refractive  dielectric  elements  in  air  as  separating  layers  between  graphene 
sheets.  Such  separating  layers  will  have  effective  medium  parameters  close  to  those  of  free  space. 
Recent  developments  in  graphene  shape  control  [14]  indicate  that  the  manufacturing  and  handling 
of  graphene  has  reached  the  stage  where  the  creation  of  the  above  stacks  is  possible.  For  simplicity 
we  treat  the  layers  between  the  graphene  sheets  as  free  space  when  calculating  the  effective  medium 
parameters.  Assuming  also  sheet  conductivity  of  uq  and  inter-sheet  distances  d  that  are  within  the 
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Figure  7:  (a)  Sketch  of  the  graphene  stack  structure,  (b)  Normalized  reflected  and  transmitted  power  for 
a  graphene  stack  illuminated  by  a  plane  wave  in  a  frequency  corresponding  to  red  light  (A  =  630  nm)  as  a 
function  of  the  stack  thickness  D  for  different  inter-sheet  distances  d. 

above  range  of  meaningful  homogenization,  we  find  that 

,  2 i  2 k0d  a0  .  . 

€t  =  1  H - ,  where  7  = -  and  a  = - .  (2) 

7  na  7T£o  c 

Here,  et  is  the  transversal  component  of  the  anisotropic  relative  effective  permittivity  and  a  is  the  fine 
structure  constant.  The  other  relative  effective  parameters  of  the  stack  are  the  normal  component 
of  the  permittivity  en  =  1  and  the  permeability  fi  =  1.  The  anisotropic  imaginary  part  of  the 
transversal  permittivity  can  be  easily  understood  as  stemming  from  the  conductivity  of  the  sheets  in 
the  transverse  direction.  Electric  fields  normal  to  the  sheets  will  not  induce  currents  and  thus  only 
the  transverse  component  of  the  permittivity  is  affected  by  the  existence  of  the  sheets. 

Figure  7(b)  shows  the  normalized  reflected  and  transmitted  power  for  a  graphene  stack  illumi¬ 
nated  by  a  plane  wave  in  a  frequency  corresponding  to  red  light  (A  =  630  nm)  as  a  function  of  the 
stack  thickness  (D)  for  stacks  of  different  inter-sheet  distances  (d).  It  is  evident  that  as  the  thickness 
of  the  stack  (D)  increases,  the  dominant  portion  of  power  that  is  not  absorbed  in  the  stack  is  due 
to  reflection  at  the  interface  of  the  stack.  Next,  we  compare  those  results  to  the  theoretical  results 
given  in  [15]  for  p-polarized  incident  plane-waves  impinging  a  medium  comprising  an  array  of  aligned 
nanotubes  of  10  nm  diameter  and  having  a  50  nm  array  spacing.  In  order  to  achieve  the  same  low 
reflectivity  (0.02%)  obtained  for  the  nanotube  array  using  the  stack,  we  need  to  place  the  sheets 
at  a  distance  of  d  =  41.2  nm.  This  yields  an  absorption  coefficient  5  times  larger  than  that  of  the 
nanotube  array.  Such  a  large  absorption  coefficient  implies  that  we  could  use  a  much  thinner  coating 
in  order  to  achieve  properties  similar  to  those  of  the  nanotube  array.  By  increasing  the  inter-sheet 
distance  (d),  the  reflectivity  is  reduced,  however,  this  has  the  detrimental  effect  of  a  decrease  in 
the  absorption  coefficient.  Thus,  by  increasing  the  inter-sheet  distance  (d),  together  with  the  entire 
thickness  of  the  stack,  one  can  increase  the  blackness  of  the  resulting  stack  without  bound.  Note 
that  when  d  is  no  longer  small  compared  to  the  wavelength,  the  homogenization  and  effective  prop¬ 
erties  we  obtained  are  no  longer  valid.  However,  we  suspect  that  the  physical  mechanism  making  the 
material  black  will  still  make  the  materials  darker,  and  we  plan  to  study  the  interpretation  when  not 
in  the  homogenization  regime.  Another  way  to  increase  the  material  blackness  would  be  to  reduce 
the  reflectivity  at  the  first  few  layers  and  absorb  the  light  in  the  bulk  of  the  stack.  To  this  end,  we 
suggest  to  continuously  reduce  the  inter-sheet  distance  (d)  with  spatial  location,  and  to  optimize  the 
functional  change  in  d  to  improve  the  blackness. 

We  are  working  with  Professor  Zhihong  Chen  at  Purdue  and  one  of  her  students  on  a  fabrication 
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method  for  the  graphene  stack.  Beyond  black  materials,  this  demonstration  will  provide  a  high-speed 
photodetector  material  and  may  also  prove  valuable  in  electrodes  for  solar  cells. 

4  Optical  and  Terahertz  Magnetic  Materials 

We  suggest  a  cellular  structure  where  each  cell  comprises  a  single  carbon  nanotube  (CNT)  coil,  and 
each  pair  of  neighboring  cells  are  of  different  handedness  (i.e. ,  left-hand  and  right-hand  coils,  as 
in  Fig.  8(a)).  This  structure  is  on  average  symmetric  with  respect  to  mirror  reflection,  and  thus 
non-chiral.  We  have  studied  the  homogenized  permeability  of  a  supercell  of  8  CNW  coils,  shown  in 
Fig.  8(b),  in  a  doubly  periodic  lattice.  We  used  a  vector  electromagnetic  solution  and  a  parameter 
extraction  procedure  similar  to  the  one  used  in  [16],  as  a  function  of  material  thickness,  in  order 
to  establish  converged  bulk  parameters.  To  simplify  our  initial  calculations  and  obtain  results  that 
are  dependent  only  on  the  ratio  between  dimensions  and  wave-length,  we  used  classical  graphite 
properties  for  the  coils  and  assumed  a  small  skin  depth  compared  with  coil  wire  thickness.  This 
simplification  allowed  us  to  concentrate  on  the  choice  of  coil  and  cell  dimensions  that  will  decrease 
the  resonant  frequency  while  keeping  the  coil  to  coil  interaction  weak,  so  that  the  effective  medium 
parameters  will  converge  even  for  small  material  thickness.  The  result  in  Fig.  8(c)  shows  the  real 
and  imaginary  components  of  the  homogenized  permeability  as  having  a  resonance,  and  that  it  is 
a  metamaterial.  It  would  therefore  appear  that  a  mixture  of  CNT  coils  could  provide  negative  /U 
at  THz  frequencies  with  an  appropriate  geometry.  We  also  found  that  the  resonant  frequency  in 
the  homogenized  result  is  insensitive  to  source  polarization  and  orientation,  which  is  very  important 
for  random  mixtures.  Our  results  on  magnetism  from  CNT  coils  are  being  developed  for  a  journal 
paper. 


(a)  (b)  (c) 

Figure  8:  (a)  Left-handed  and  right-handed  carbon  coils  that  are  placed  in  a  x  a  x  a  cells.  The  coil  diameter 
and  length  is  0.5a,  and  the  wire  diameter  is  0.02a.  (b)  Non-chiral  unit-cell  of  periodic  metamaterial  that 
comprises  4  left-handed  coils  and  4  right-handed  coils,  (c)  Real  (//)  and  imaginary  (//")  components  of  the 
effective  permeability  of  the  coil  metamaterial,  obtained  using  a  scattering  model. 


5  Electromagnetic  Energy  in  Dispersive  Materials 

A  negative  refractive  index  can  be  achieved  when  both  the  dielectric  constant  e  and  the  permeability 
fj,  are  negative,  which  requires  an  operating  frequency  near  a  resonance  in  both  constitutive  param¬ 
eters.  In  this  frequency  regime,  the  material  properties  can  be  highly  dispersive  for  the  candidate 
metamaterials,  and  the  loss  can  be  substantial.  Loudon  [17],  Ruppin  [18],  Tretyakov  [19],  Boardman 
and  Marinov  [20],  Cui  and  Kong  [21],  and  Nunes  et  al.  [22]  and  several  others  have  studied  energy 
density  issues  in  dispersive  media.  This  substantial  body  of  work  presents  results  where  all  energies 


(total,  stored  and  dissipated)  remain  positive.  However,  Ziolkowski  [23]  questions  whether  the  elec¬ 
tromagnetic  energy  must  remain  positive.  We  have  shown  that  the  electric  field  energy  density  in 
dispersive  materials  can  be  exactly  separated  into  stored  and  lost  energies  when  using  the  electric 
field  as  the  basis  [24],  Our  analytical  and  numerical  results  for  a  narrowband  modulated  wave  near  a 
material  resonance  indicated  that  the  stored  energy  can  be  negative,  while  the  total  energy  remains 
positive  in  this  high-loss  regime.  We  have  obtained  the  energy  density  for  a  negative  refractive  index 
material.  For  the  special  case  of  a  frequency-independent  wave  impedance  in  the  material,  the  exact 
separation  of  both  electric  and  magnetic  energy  densities  into  stored  (electric  and  magnetic)  and  lost 
(electric  and  magnetic)  densities  is  shown  to  be  possible. 

Using  Poynting’s  theorem,  the  time  derivative  of  electric,  ue ,  and  magnetic,  uh,  field  energy 
densities  can  be  expressed  as 


8ue 

~df 


=  E  •  — —  and 
dt 


duH  _  H  dB 
dt  ~  dt' 


(3) 


Expressing  the  electric  flux  density  as  an  exact  Fourier  superposition  using  the  complex  material 
constitutive  parameters  and  choosing  the  electric  field  as  the  basis  (see  [24]  for  details),  we  have 
shown  that  the  electric  field  energy  density  can  be  separated  into  the  stored  and  lost  electric  field 
energy  densities.  Similarly  for  the  magnetic  field  case,  with  magnetic  the  field  as  the  basis,  the 
magnetic  field  energy  density  has  been  separated  into  the  stored  and  lost  magnetic  field  energy 
densities  [25]. 

The  electric  and  magnetic  fields  are  related  through  H{t)  =  (l/27r)  f/^oE((v)r/~1(uj)e~iU,tduj, 
where  iq{uj)  is  the  frequency-dependent  wave  impedance.  When  the  electric  field  is  chosen  as  the 
basis,  ue  is  exactly  separable  into  we  (stored  electric  field  energy)  and  qE  (lost  electric  field  energy), 
but  uh  is  generally  not.  Similarly,  when  the  magnetic  field  is  chosen  as  the  basis,  uh  is  exactly 
separable.  However,  in  the  special  case  when  both  e(u)  and  n(ui)  are  constant  or  e(uj)  =  /u(w), 
rj(io)  is  independent  of  frequency,  giving  H(t )  =  r/E(t),  where  rj  is  the  constant  value  of  the  wave 
impedance  in  the  frequency  domain.  For  this  case,  the  total  energy  density,  u  =  ue  +  uh,  is  exactly 
separable  into  stored  (which  is  returned  to  the  field,  w  =  we  +  wh)  and  lost  (converted  to  a  non- 
electromagnetic  form,  q  =  qE  +  QH )  energies. 

We  assume  E  =  e(t)cos(uiot) ,  with  slowly  varying  modulation  signal  e(t)  relative  to  to  =  ‘Ztt/ojq, 
producing  an  effective  bandwidth  that  is  small  relative  to  the  features  of  e(uj).  Following  our  earlier 
work  [24],  the  approximate  form  of  dD/dt  is  obtained  by  using  the  first  two  terms  in  the  Taylor 
series  expansion  of  (u>  +  uJo)e(co  +  u>o)  about  u  =  (1.  The  product  of  E(t)  and  the  approximate  form 
of  dD/dt  gives  the  approximate  form  of  duE/dt.  After  some  additional  steps  (see  [24]),  we  obtain 
the  average  electric  field  energies 


{we) 


( 


dqE 

dt 


1  d  [uie'{uj)\ 
4£°  du 

■^uj0e0e"e2(t). 


e\t) 


(4) 

(5) 


Following  a  similar  analysis  using  the  magnetic  field  as  the  basis  instead  of  the  electric  field,  we 
obtain  the  average  magnetic  field  energies 
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where  h(t)  is  the  slowly  varying  envelope  of  the  magnetic  field.  Equations  (4)  -  (7)  apply  regardless 
of  the  degree  of  loss,  provided  that  the  two-term  Taylor  series  expansion  for  dD(t)/dt  and  dB(t)/dt 
are  sufficiently  accurate. 

For  the  simple  case  when  rj  is  independent  of  frequency  and  the  energy  density  is  exactly  sepa¬ 
rable,  we  choose  identical  electric  and  magnetic  susceptibilities  ( xe  =  Xm)  in  arriving  at  a  negative 
index  material  given  by  Xe(w)  =  ci\  (u>2  —  uo2  —  i'yiu)  1 .  Example  dispersive  dielectric  constant  and 
relative  permeability  responses  are  shown  in  Fig.  9(a).  The  material  parameters  are  chosen  so  that 
n'iujo)  =  —  1  for  u>o  =  8/7,  where  n'{ ojq)  is  the  real  part  of  the  refractive  index  of  the  material. 

We  assume  a  Gaussian  modulated  field  given  by  E(t)  =  (oV 2tt)~  1  exp  [~(t  —  tc)2(2a2)  x]  cos[cuo(i— 
fc)]  with  lo o  =  4/9,  8/7  and  4/3,  tc  =  210(27t),  a  =  27(27t),  and  carrier  phase  delay  wq tc.  The  magnetic 
field  was  obtained  from  the  Fourier  decomposition  of  this  electric  field  using  H(uj)  =  E(u)/ri(u). 
The  expression  for  duE/dt  in  (3)  was  evaluated  numerically,  with  D(t)  formed  by  Fourier  trans¬ 
forming  (using  an  FFT)  the  temporal  electric  field,  multiplying  by  e(cu),  and  then  inverse  Fourier 
transforming  the  result.  The  quantities  dwE/dt  and  dqE/dt  were  formed  using  the  exact  decom¬ 
positions.  Numerical  integration  over  time  followed  by  time-averaging  over  the  period  to  resulted 
in  (we),  ( dqE/dt )  and  (ue)-  Using  the  magnetic  field  instead  of  the  electric  held  and  following  a 
similar  procedure,  the  average  quantities  ( wh ),  ( dqu/dt )  and  (uh)  were  obtained. 

Figure  9(b)  shows  that  when  the  carrier  frequency  is  sufficiently  far  from  the  resonant  frequency, 
(we)  as  well  as  (wh)  are  positive,  just  as  they  would  have  been  in  free  space.  However,  when  the 
carrier  frequency  is  just  above  the  resonant  frequency  (cuo  =  8/7),  and  in  this  case  where  the  refractive 
index  is  negative,  (we)  as  well  as  (wh)  become  negative.  Figure  9(c)  shows  the  numerical  results 
for  ( dqE/dt )  and  (dqu/dt)  obtained  for  the  same  three  carrier  frequencies,  along  with  the  results 
from  the  approximate  analytical  expressions  in  (5)  and  (7).  Throughout,  the  total  average  energy 
(( ue)  +  (uh ))  remains  positive,  as  Fig.  9(d)  shows.  At  times  long  after  the  pulse  has  gone,  (qs)  =  (ue) 
and  (qu)  =  (uh),  indicating  that  no  additional  energy  beyond  that  delivered  by  the  electromagnetic 
held  is  involved,  despite  (we)  and  (wh)  being  negative  for  some  time.  Notice  the  excellent  agreement 
between  the  exact  results  and  the  approximate  analytical  expressions  in  Figs.  9(b)  and  9(c). 

The  energies  represent  running  totals  as  a  function  of  time,  and  the  sign  of  dw/dt  gives  the 
transfer  direction  of  the  transient  energy,  which  integrated  over  time  and  averaged  over  the  carrier 
period  to  give  the  net  energy  exchange  up  to  each  instant  in  time.  When  (dw/dt)  is  positive  during 
the  first  half  of  the  Gaussian  and  negative  during  the  second  half,  (w)  in  Fig.  9(b)  is  positive 
throughout  and  zero  for  large  time,  after  the  pulse  has  left  that  point  in  space.  In  the  case  where 
we  find  negative  stored  energy,  near  a  resonance,  (dw/dt)  is  first  negative,  and  then  positive.  The 
negative  field  energy  is  energy  borrowed  from  the  material,  and  during  the  second  half  of  the  pulse, 
the  field  repays  this  energy. 

6  Electromagnetic  Plane  Wave  Forces  in  Homogeneous  Materials 

We  consider  the  various  contributions  to  the  electromagnetic  force  from  a  fundamental  perspective, 
including  the  material  constitutive  parameters  and  their  dispersive  properties,  and  show  that  the 
radiation  pressure  is  positive  for  propagating  waves  in  homogeneous  passive  media,  even  when  the 
refractive  index  is  negative  [26,27].  We  show,  however,  that  the  force  can  be  negative  when  the 
medium  has  gain,  and  we  describe  an  experiment.  Our  work  provides  a  foundation  for  the  treatment 
of  forces  in  dispersive  inhomogeneous  material  systems,  such  as  objects  in  some  background  material, 
and  the  further  study  of  the  influence  of  the  background  material  properties. 
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Figure  9:  (a)  Real  (solid  line)  and  imaginary  (dotted  line)  part  of  the  dielectric  constant  and  rel¬ 
ative  permeability  {e,  fi}  =  1  +  —  lo2  —  iO.Su;)-1  for  uj\  =  1.  The  three  different  carrier 

frequencies  used  (u)q  =  4/9, 8/7, 4/3)  are  marked  by  thin  vertical  solid  lines,  (b)  (ue)  \e"=o  =  (we) 
and  {uh)\h"=o  =  (wh)  obtained  for  {e',//}  =  Real{e, /r}.  The  lines  give  exact  numerical  results 
and  the  model  result  (circles)  plots  (4)  for  the  luq  =  8/7  case,  (c)  (Oue /dt)\e'=o  =  (dqE/dt)  and 
(dun /dt)\fj,i=o  =  ( dqu/dt )  obtained  with  =  Imag{e,  ^}.  The  lines  give  exact  numerical 

results  and  the  model  result  (circles)  is  from  (5)  for  the  luq  =  8/7  case,  (d)  Exact  time-averaged 
energy  densities,  (ue)  and  (uh),  obtained  for  {e,  /r}. 
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The  electromagnetic  plane  wave  force  density  in  source-free  homogeneous  medium  is  [26] 


f 

For  a  modulated  electric  field  E(f) 
force  density  to  be 


op  dfi0M 

at"  ,loH  at~  e“E' 


(8) 


=  e  E(t)  =  e  e(t)  cos(ujq t),  we  have  shown  the  normalized  average 


(fn) 


(C\  +  D\)  e2(t )  +  ( C2  +  D2)  e(t) 


d  e(t) 
dt 


+  C3 


<9e(f)l  2 
dt 


(9) 


where  the  constants  C’s  and  D’s  depend  on  the  constitutive  material  parameters  and  are  defined  in 
Ref.  [26]. 

To  evaluate  the  significance  of  the  force  in  a  homogeneous  medium  with  loss  and  gain  medium, 
consider  a  Gaussian  modulation  signal  given  by  e(t)  =  exp[~ t2  (2cr2)  ],  where  we  use  a  =  287r, 

ojq  =  1,  and  N  =  220  sample  points  over  a  temporal  support  of  2137t.  A  numerical  simulation 
was  carried  out  by  first  evaluating  the  electric  field,  then  determining  the  magnetic  field  using 
a  discrete  Fourier  transform  and  by  incorporating  the  frequency-dependent  material  parameters. 
Subsequently,  the  total  instantaneous  plane  wave  force  was  obtained  using  the  scalar  form  of  (8),  and 
this  was  numerically  integrated  over  time  to  form  (/).  We  describe  this  evaluation  as  exact  because 
the  accuracy  is  subject  only  to  numerical  precision.  The  time-average  force  was  also  calculated 
using  the  analytical  model  given  in  (9).  To  evaluate  the  coefficients  C2,  C3,  and  D2 ,  we  used 
numerical  derivatives  of  fi(u)  and  e(u;).  Figures  10(a)  and  (b)  show  the  time-averaged  normalized 
electromagnetic  plane  wave  force  density  calculated  using  the  exact  numerical  procedure  (solid  line) 
and  the  analytical  expression  in  (9)  (circles)  for  a  non-magnetic  material  with  /i  =  1  and  having 
overall  (a)  loss,  with  permittivity  eioss,  and  (b)  gain,  with  permittivity  egain.  We  set  eioss  =  e  = 
1  +  2(0. 952  —  u)2  —  i0.2 for  a  passive  and  lossy  two-level  system  and  egain  =  e  =  1  —  0.6(1.052  — 
cj2  —  i0.2o;)_1  for  an  active  system  having  overall  gain.  From  Fig.  10,  (fn)  is  positive  for  the  lossy 
material  and  negative  for  the  material  with  gain,  throughout  the  pulse.  Notice  that  the  analytical 
result  of  (9)  is  in  excellent  agreement  with  the  numerical  data. 

In  a  physical  demonstration  of  a  negative  force,  the  gain  medium  would  be  excited  to  create 
a  population  inversion,  and  this  system  would  be  illuminated  with  a  probe  electromagnetic  wave 
[26,27]..  The  force  will  be  in  a  direction  opposite  to  the  Poynting  vector  for  this  probe  wave.  To 
evaluate  the  significance  of  the  negative  plane  wave  force  in  media  with  gain,  we  compared  the 
expected  positive  force  for  a  metal,  which  has  been  measured  [28],  with  the  estimated  negative 
force  in  an  excited  quantum  dot  medium  [4],  A  negative  force  on  a  material  exhibiting  gain  can  be 
considered  a  collective  result  of  the  stimulated  emission  from  each  atom,  where  the  atomic  recoil  is 
in  the  direction  opposite  to  that  of  the  emitted  photons  [29] . 


7  Anisotropic  Metamaterials  Based  on  Metal-Insulator  Stacks 

Multi-layer  metal-insulator  stacks  at  infra-red  and  visible  frequencies  can  be  characterized  as  uniax- 
ially  anisotropic  metamaterials,  provided  the  periodicity  is  much  smaller  than  the  operating  wave¬ 
length.  Effective  medium  parameters  characterizing  the  anisotropic  media  typically  involve  two 
variables:  ej_,  the  permittivity  component  perpendicular  to  the  material  surface,  and  en,  the  com¬ 
ponent  parallel  to  the  surface.  It  has  been  shown  that  the  strong  anisotropy  can  lead  to  a  wide 
variety  of  physical  phenomena  and  important  applications.  Here  we  describe  two  relevant  results: 
subwavelength  imaging  and  a  miniaturized  spectrometer. 
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Figure  10:  The  time-averaged  normalized  electromagnetic  plane  wave  force  density  calculated  using 
the  exact  numerical  procedure,  shown  by  solid  lines,  and  the  analytical  expression  in  (9),  shown  by 
circles,  for  a  non-magnetic  material  with  /U  =  1  and  having  overall  (a)  loss,  with  permittivity  eioss 
and  (b)  gain,  with  permittivity  egain.  The  electric  field  modulation  is  given  by  eg(t),  and  the  carrier 
frequency  is  ujq  =  1. 


We  have  numerically  shown  a  bilayer  lens  with  subwavelength  resolution  which  consists  of  one 
anisotropic  layer  and  one  isotropic  dielectric  layer,  with  leaky  waveguide  operation  of  the  anisotropic 
film  [30].  In  Fig.  11(a),  we  show  the  homogenized  Ag/Si02  metal-insulator  stack  on  a  silicon  (Si) 
film  as  an  example.  The  finite  element  simulation  result  clearly  demonstrates  that  the  slab  behaves 
as  an  optical  leaky  waveguide  radiating  backwards.  Because  the  leaky  mode  supported  by  the 
system  can  have  |Re(/cz)|  >  ko,  a  planar  near-field  imaging  system  that  transfers  sub-diffraction 
features  can  be  envisioned  by  pairing  a  backward  leaky  waveguide  with  a  silicon  slab,  as  shown  in 
Figs.  11(b)  and  11(c).  As  illustrated  in  Fig.  11(b),  the  waves  radiated  from  the  source  emerge  as 
strongly  confined  beams  inside  the  anisotropic  slab,  and  then  re-radiate  into  the  backward  quadrant 
at  the  interface  between  the  anisotropic  waveguide  and  the  silicon  slab  to  form  a  focal  point  at  a 
certain  distance  from  the  interface.  In  contrast  to  the  realization  of  a  super-resolving  lens  using 
stack-based  anisotropic  bilayer  lenses  that  we  have  analyzed  and  developed  previously  [31],  the 
present  approach  to  subwavelength  imaging  does  not  impose  any  restrictions  on  impedance  match 
or  operating  wavelength,  and  it  requires  the  synthesis  of  only  one  anisotropic  bulk  medium. 

We  have  designed  a  new  spectrometer  class  that  could  be  realized  with  TM  illumination  of  a  small 
aperture  placed  in  front  of  a  metal-insulator  stack,  and  we  have  derived  analytical  formulations  to 
evaluate  the  working  wavelength  range  and  the  wavelength  sensitivity  [32].  This  spatial  filtering 
property  is  based  on  the  frequency-dependent  preferred  propagation  direction  associated  with  an 
anisotropic  slab,  determined  by  — e^/e^,  due  to  intrinsic  material  dispersion.  As  an  example,  we 
investigated  a  Ag/ZnO  metal-insulator  stack  as  a  possible  realization  of  exez  <  0  for  the  spectrometer. 
Figure  12(a)  depicts  the  frequency  response  of  the  homogenized  Ag/ZnO  stack  at  the  output  plane  for 
three  different  wavelengths,  and  compares  the  magnitude  of  the  transmitted  electric  fields  employing 
the  closed- form  solution  we  developed  (for  design  purposes),  which  we  derived  based  on  a  Green’s 
function  analysis  [32],  and  full- wave  finite  element  simulation  results.  In  Fig.  12(b),  we  plot  the 
position  of  the  major  peak  as  a  function  of  wavelength,  showing  that  the  principal  peak  continuously 
shifts  away  from  the  center  as  wavelength  increases,  and  that  a  larger  slab  thickness  gives  rise  to 
larger  peak  shifts. 
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Figure  11:  (a)  i/-field  distribution  of  a  Ag/SiC>2  bulk  anisotropic  leaky  waveguide  embedded  in  a  Si  host. 
The  volume  fraction  for  Ag  is  0.2.  The  width  of  the  slab  is  0.16 /im.  A  magnetic  current  line  source  is  located 
at  (x  =  0,  z  =  —  0.25/rm).  A  =  563. 6nm,  esi  =  16.3  +  *0.26,  £Ag  =  — 11.89  +  *0. 83,  esio2  =  2.13,  giving 
ex  =  —0.67  +  *0.17,  ez  =  2.79  +  *9.05  x  10-3.  (b)  A  Ag/Si02  bulk  anisotropic  leaky  waveguide  juxtaposed 
with  a  Si  slab,  functioning  as  a  sub-diffraction  imaging  system.  The  thickness  of  the  waveguide  d  =  0.16  fim, 
the  thickness  of  silicon  slab  L  =  0.45  fim,  and  the  surrounding  medium  is  free  space.  The  excitation  current 
line  source  is  located  at  x  =  0.  (c)  The  field  profile  at  the  image  plane  x  =  d  +  L. 
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Figure  12:  An  anisotropic  spectrometer,  (a)  Electric  field  magnitude  at  three  wavelengths  obtained 
from  finite-element  simulation  and  analytical  results:  a  magnetic  line  current  source  radiates  through 
an  anisotropic  Ag/ZnO  slab  with  D  =  0.5,  and  the  slab  thickness  d  =  100  nrn.  (b)  Variation  of  the 
position  the  principal  peak  as  a  function  of  wavelength  for  a  homogenized  Ag/ZnO  stack  slab. 


We  have  investigated  wave  propagation  behavior  inside  a  cylindrically  anisotropic  metamate¬ 
rial  with  a  general  parameter  space  <  0,  and  identified  the  shape  of  the  resonance  cones  as 
two  sets  of  logarithmic  spiral  curves  [33].  The  formation  of  resonance  cones  inside  a  curved  film 
of  cylindrically  anisotropic  material  excited  by  a  line  source  is  applied  to  demonstrate  a  new  spec¬ 
trometer/multiplexer  class  that  exploits  metal-insulator  stack  metamaterials.  A  bilayer  imaging 
system  consisting  of  two  cylindrically  curved  shells  of  such  media,  which  holds  promise  for  far-held 
sub-diffraction-limit  imaging  and  obviates  the  strict  condition  of  material  parameters,  has  also  been 
proposed  [34]. 
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8  Impact  of  Surface  Roughness  on  the  Effective  Medium  Parame¬ 
ters  and  Subwavelength  Imaging 


The  effective  dielectric  constant  for  a  multilayer  metal-insulator  stack  is  obtained  from  a  numeri¬ 
cal  simulation  and  compared  with  the  analytical  result.  The  impact  of  the  surface  roughness  for 
subwavelength  imaging  is  established  [35].  This  work  is  now  being  prepared  as  a  journal  paper. 

The  evolution  of  nanofabrication  techniques  has  allowed  the  fabrication  of  precise  structures. 
However,  the  metal-insulator  stacks  that  can  yield  far-subwavelength  image  resolution  (see,  for  ex¬ 
ample,  [1,36])  require  very  smooth  and  very  thin  metal  and  insulator  films.  It  is  of  importance  to 
understand  the  influence  of  this  inevitable  roughness  for  the  improved  design  of  nanostructures  and 
experiments.  Here,  we  use  a  simple  method  for  the  effective  medium  parameter  extraction  for  a 
metal-insulator  planar  slab  and  study  the  effect  on  the  subwavelength  imaging  performance  using 
silver  (Ag)  and  zinc  oxide  (ZnO)  as  representative  materials. 

The  numerical  simulation  used  for  the  extraction  procedure  is  carried  out  using  an  FEM  solver 
(COMSOL  Multiphysics  [37]).  The  geometry  consists  of  a  100  nm  thick  planar  slab  placed  in  vacuum 
and  composed  of  alternate  layers  of  Ag  and  ZnO.  The  two  ports  (left-  and  right-most  boundaries) 
are  located  at  a  distance  of  600  nm  from  each  slab  surface.  The  top  and  bottom  boundaries  are 
perfect  electric  conductors  with  a  separation  of  1  //rri,  forming  a  waveguide,  and  a  TM  wave  is 
incident  normally  on  the  slab  from  the  left  side.  The  Ag  and  ZnO  dielectric  constants  {(-Ag  and  eznOi 
respectively)  are  obtained  by  interpolating  measured  bulk  data  [38,39].  The  scattering  parameters  at 
the  ports  are  obtained  from  the  FEM  simulation,  and  using  these,  the  effective  transverse  dielectric 
constant  is  calculated  according  as  [13] 


eeff  = 


(1  -  Sii)2  -  S, 


21 


(l  +  Sn)2-S2V 


Mf  ~  L 


(10) 


We  have  assumed  \ieff  =  1  f°r  °ur  numerical  simulations  because  the  materials  are  non-magnetic.  In 
the  homogenization  limit,  the  effective  transverse  dielectric  constant  for  a  planar  slab  with  smooth 
surfaces  depends  on  the  duty  cycle,  D ,  and  is  analytically  given  by  [40] 


(eff  =  £AaD  +  ezno{^  —  D).  (11) 

Figure  13(a)  shows  eeff  obtained  numerically  for  a  planar  100  nm  thick  Ag/ZnO  multilayer 
slab  for  different  numbers  of  periods,  and  compared  with  the  analytical  homogenized  result  in  (11). 
It  is  observed  that,  in  this  case,  10  periods  is  sufficient  to  approach  the  homogenization  limit.  Fig¬ 
ures  13(b)  and  (c)  show  the  effect  of  the  surface  roughness  on  eeff  for  the  100  nm  Ag/ZnO  multilayer 
slab  having  5  and  10  periods,  respectively.  As  we  increase  the  surface  roughness,  the  deviation  of 
eeff  extracted  from  the  numerical  model  increases,  relative  to  that  determined  by  the  analytic  result 
in  (11).  The  standard  deviation  in  eeff  also  increases  with  increase  in  surface  roughness.  It  is  thus 
clear  that  roughness  has  a  substantial  influence  on  the  homogenized  dielectric  constant  determined 
using  the  FEM  model,  which  provides  data  that  could  be  representative  of  a  measurement. 

To  evaluate  the  impact  of  surface  roughness  in  an  imaging  application,  we  perform  another  series 
of  numerical  studies.  Figure  14(a)  shows  a  stack  having  statistical  roughness,  and  from  a  series  of 
such  calculations,  we  generated  Fig.  14(b),  which  shows  the  impact  on  subwavelength  resolution. 
These  calculations  assumed  a  standard  deviation  of  around  1  nm,  comparable  with  that  achieved 
with  e-beam  evaporation  and  a  35  nm  him  thickness  [41]. 

This  work  constitutes  another  step  towards  a  better  design  of  manufacturable  stack-based  nanos¬ 
tructures  by  improving  our  understanding  of  the  impact  of  the  surface  roughness  on  the  effective 
dielectric  constant.  The  metal-insulator  stack  shows  promise  for  subwavelength  imaging  and  other 
applications. 


15 
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Figure  13:  The  real  (top  row)  and  imaginary  (bottom  row)  part  of  the  effective  transverse  dielectric  constant,  ee// 
obtained  as  a  function  of  wavelength  for  the  100  nm  Ag/ZnO  multilayer  slab,  (a)  ee//  for  different  numbers  of  periods 
with  smooth  surfaces,  (b)  and  (c)  ee//  for  the  slab  with  5  and  10  periods,  respectively,  obtained  for  different  standard 
deviations  of  the  surface  roughness,  a.  The  lines  and  the  error  bars  show  the  mean  and  the  standard  deviation  for  the 
extracted  ee//-  The  analytical  result  shown  in  green  is  obtained  using  (11),  while  the  others  are  obtained  from  the 
FEM  simulation  and  using  (10). 
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Figure  14:  (a)  The  geometry  considered  for  our  numerical  study  of  the  influence  of  surface  roughness 
in  an  anisotropic  lens:  a  screen  with  two  10  nrn  slits  and  30  nm  center-to-center  distance;  Ag/ZnO 
multilayer  stack  with  5nm  thick  layers  and  a  surface  roughness  with  standard  deviation  cr;  and 
incident  field  wavelength  A  =  413.3  nm.  (b)  Statistical  influence  of  random  film  roughness  on  the 
resolution  in  a  double  slit  imaging  arrangement,  showing  that  achievable  film  smoothness  should 
allow  far  subwavelength  image  resolution. 
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